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We apply neutron diffraction, high-resolution synchrotron x-ray diffraction, magnetization mea- 
surements, electronic structure calculations, and quantum Monte-Carlo simulations to unravel 
the structure and magnetism of (CuCl)LaTa207. Despite the pseudo-tetragonal crystallographic 
unit cell, this compound features an orthorhombic superstructure, similar to the Nb-containing 
(CuX)LaNb207 with X = CI and Br. The spin lattice entails dimers formed by the antiferromag- 
netic fourth-neighbor coupling J4, as well as a large number of nonequivalent interdimer couplings 
quantified by an effective exchange parameter J e g. In (CuCl)LaTa207, the interdimer couplings are 
sufficiently strong to induce the long-range magnetic order with the Neel temperature Tjv — 7 K 
and the ordered magnetic moment of 0.53 [is, as measured with neutron diffraction. This magnetic 
behavior can be accounted for by J e ff/J4 — 1-6 and J4 ^ 16 K. We further propose a general mag- 
netic phase diagram for the (CuCl)LaNb2 07-type compounds, and explain the transition from the 
gapped spin-singlet (dimer) ground state in (CuCl)LaNb207 to the long-range antiferromagnetic 
order in (CuCl)LaTa207 and (CuBr)LaNb207 by an increase in the magnitude of the interdimer 
couplings J e fi/J4, with the (CuCl)LaM2 07 (M = Nb, Ta) compounds lying on different sides of the 
quantum critical point that separates the singlet and long-range-ordered magnetic states. 



PACS numbers: 61.66.Fn, 75.30.Cr, 75.30.Et, 75.10.Jm 

I. INTRODUCTION 

Quantum phase transitions are one of the hot topics 
in present-day solid-state physics. Experimental stud- 
ies in this field enable a test of existing theories on 
quantum criticality as well as a direct access to exotic 
phases emerging in the vicinity of quantum critical points 
(QCPs). This task, however, remains highly challeng- 
ing, because model systems amenable to the experimen- 
tal study should lie close to the QCP and allow for a 
continuous tuning across the QCP. Such tuning can be 
achieved by applying external pressure or performing 
chemical substitutions. While external pressure facili- 
tates continuous evolution of the system, the pressure 
cell complicates experimental work and restricts the set 
of applicable experimental techniques. Although mea- 
surements on chemically substituted samples are experi- 
mentally more feasible, the downside of the chemical ap- 
proach is the inevitable disorder that may be detrimental 
for the physical effect under consideration. 

Magnetic systems have become one of the main play- 
grounds for experimental and theoretical studies of quan- 
tum critical phenomenal— In insulating magnets, the ef- 
fect of chemical substitution strongly depends on the po- 
sition of the replaced atom and its contribution to mag- 
netic couplings. Replacement of a ligand (anion) will 
usually disrupt the spin lattice, because ligand orbitals 
are an integral part of superexchange pathways that are 
responsible for the coupling. While systems with dis- 
rupted spin lattices are quite interesting on their owe^- 
effects of bond randomness are highly unfavorable for ex- 
perimental studies of quantum phase transitions in per- 
fect, non-disrupted spin lattices. Cation replacement is 



a better method for tuning the system, because the spin 
lattice remains nearly homogeneous, whereas individual 
exchange couplings arc only slightly modified owing to 
the change in lattice parameters, electrostatic fields, and 
crystal-field splittings. However, the effect of the cation 
replacement is often too weak to modify the magnetic 
ground state and drive the system across the QCP. 

The (CuCl)LaNb2- x Ta x 07 solid solutions are one of 
the promising systems with a possible quantum phase 
transition induced by the cation replacement^ The Nb 
compound [x = 0) reveals a singlet ground state with a 
spin gap^ whereas the Ta compound (x = 1) is long- 
range antiferromagnctically ordered below T/v — 7 K. 11 
The compositions with fractional x values are intermedi- 
ate between these two different regimes and presumably 
form a combination of the gapped singlet and gaplcss 
long-range-ordered magnetic phases 

Crystal structures of the (CuX)LaM 2 07 compounds 
(X = CI and Br, M = Nb and Ta) feature mag- 
netic [CuX] layers sandwiched between non-magnetic 
[LaM2C>7] perovskite slabs . 13 ' 14 Previous studies of the 
Nb-based systems established the formation of an or- 
thorhombic superstructure related to the Jahn- Teller dis- 
tortion of Cu +2 , the ordering of X atoms, and the tilt- 
ing distortion within the [LaNb2C>7] slabs* 1 ^— The re- 
sulting spin lattice entails spin dimers with a complex 
and still controversial pattern of interdimer couplings 
in (CuCl)LaNb2 07 ) 16 ' 19 as well as strong interdimer 
couplings that trigger the long-range antiferromagnetic 
(AFM) order in (CuBr)LaNb 2 7 iA£ Surprisingly, little is 
known about the crystal structure of (CuCl)LaTa 2 07, 
and only its disordered version has been reported^ 

Because details of the low-temperature crystal struc- 
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ture and the nexus of the Ta and Nb compounds are 
crucial to understand the magnetic behavior of the 
(CuCl)LaNb2_ :c Ta 2: 07 solid solutions, we performed an 
extensive structural study of (CuCl)LaTa207. The crys- 
tallographic results were further used for a microscopic 
analysis that demonstrated a close similarity between 
the Nb and Ta systems. Combining the crystallographic 
and microscopic information, we elucidate the differ- 
ences in the atomic positions and the ensuing variation 
of the exchange couplings in the (CuX)LaM 2 Oy family. 
We provide a consistent, generalized, and quantitative 
description of the respective compounds, and address 
unresolved issues, such as the nature of magnetism of 
(CuCl)LaNb 2 7 . 



II. METHODS 

Powder samples of (CuCl)LaTa207 were prepared by 
a two-step procedure according to Ref. [HI. First, the 
RbLaTa207 precursor was obtained by firing a mixture 
of La 2 3 , Ta 2 5 , and 40 % excess of Rb 2 C0 3 at 800 °C 
for 12 h and 1050 °C for 12 h. The fused sample was 
washed with water, dried at 70 °C, and mixed with twice 
the equimolar amount of anhydrous CuCl 2 under argon 
atmosphere. This mixture was pressed into pellets, sealed 
in an evacuated quartz tube, and heated at 400 °C for 
36 h. The sample was again washed with water to remove 
the excess CuCl 2 and the RbCl formed during the reac- 
tion, and dried overnight at 70 °C. Sample purity of the 
starting materials, the RbLaTa 2 07 precursor, and the fi- 
nal product was controlled by powder x-ray diffraction 
(XRD). 

Room-temperature XRD patterns were collected with 
the laboratory Huber G670 Guinier camera (CuK Q i ra- 
diation, 28 = 3—100° angle range, ImagePlate detector). 
High-resolution XRD data were further obtained at the 
ID31 beamline of European Synchrotron Radiation facil- 
ity using a constant wavelength of A ~ 0.4 A and eight 
scintillation detectors, each preceded by a Si (111) ana- 
lyzer crystal, in the angular range 26 = 1 — 40 dcg. The 
powder sample was contained in a thin-walled borosil- 
icate glass capillary that was spun during the experi- 
ment. The sample temperature was controlled by a He- 
flow cryostat (at 4.2 K), a liquid- nitrogen cryostream 
(80 - 310 K), and a hot-air blower (350 - 750 K). 

Neutron diffraction data were collected at the high- 
resolution diffractometer D2B (A ~ 1.595 A, Q = 0.9 - 
7.6 A" 1 , T = 10 K) and the high- flux diffractometer 
D20 (A ~ 2.417 A, Q = 0.35 - 4.9 A" 1 , T = 1.5 K 
and 40 K), both installed at the Institute Laue-Langevin 
(ILL, Grenoble, France). A 6 g powder sample used for 
all neutron measurements was loaded into a cylindrical 
vanadium container and cooled down with the standard 
Orange He-flow cryostat. Crystal and magnetic struc- 
tures of (CuCl)LaTa 2 7 were refined with JANA2006 20 
and FullProf— programs, respectively. 

Thermogravimetric analysis (TGA) was performed 



with STA409 Netzsch thermal balance in the 300-1070 K 
temperature range in Ar atmosphere. The powder sam- 
ple of (CuCl)LaTa 2 07 was placed into a corundum cru- 
cible. 

The magnetic susceptibility was measured with the 
Quantum Design MPMS in the 4 — 380 K temperature 
range in the applied magnetic field of 0.5 T. 

To investigate electronic structure of (CuCl)LaTa 2 07 
and evaluate individual exchange couplings, we per- 
formed scalar-relativistic band structure calculations in 
the framework of density functional theory (DFT) with 
the local density approximation (LDA)22 and general- 
ized gradient approximation (GGA)22 for the exchange- 
correlation potential. The full-potential FPL09. 01-35 
code was used^ Reciprocal space was sampled with a 
k mesh of 256 points for the 48-atom crystallographic 
unit cell and 32 or 48 points for doubled 96-atom super- 
cells. The convergence with respect to the k mesh was 
carefully checked. The mean-field DFT+C/ correction for 
strong correlation effects was applied. Further details of 
the computational procedure are described in Sec. IIV Al 

The magnetic model derived from the DFT calcu- 
lations was further studied by quantum Monte Carlo 
(QMC) simulations using loop2£ and dirloop_sse2£ al- 
gorithms of the ALPS simulation packaged Simulations 
were performed for finite lattices with periodic boundary 
conditions. The shapes and sizes of these lattices were 
adjusted to achieve an appropriate finite-size scaling for 
the staggered magnetization (m s ) and the Neel temper- 
ature (Tjv), or to ensure the lack of finite-size effects for 
thermodynamic properties. Details of the simulations 
procedure arc described in Sections IIV CI and [V] 



III. CRYSTAL STRUCTURE 

The crystal structures of (CuCl)LaNb 2 07 and 
(CuBr)LaNb 2 07 are pseudotetragonal, with a tiny dif- 
ference between the a and b lattice parameters Al- 
though in the Nb compounds the orthorhombic splitting 
could be observed in high-resolution synchrotron XRD 
experiments, we failed to detect any signatures of a sim- 
ilar splitting in (CuCl)LaTa 2 07. The structure refine- 
ment reported below (Table [J) suggests that the differ- 
ence between the a and b parameters is indeed diminu- 
tively small and masked by the sizable reflection broad- 
ening, which is rather similar for different hkl indices and 
signifies a sizable amount of defects and/or low particle 
size. 

Owing to the large difference in the x-ray scatter- 
ing from heavy Ta and light O atoms, neutron diffrac- 
tion data were used for the structure refinement. Low- 
temperature neutron patterns (Fig. Q]) revealed 10 su- 
perstructure reflections violating the simple tetragonal 
unit cell with a su b — 3.88 A and c su b — 11.74 A, 
as reported in Ref. [l4|. Therefore, a four times larger 
2a S ub x 2a su b x c su b unit cell was used. In contrast to 
(CuBr)LaNb 2 07ji£ both superstructure reflections with 
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FIG. 1. (Color online) Low-temperature neutron diffraction 
patterns of (CuCl)LaTa207 measured at 2 K and 20 K with 
the high-flux instrument D20. Black ticks indicate the reflec- 
tions of the o su b x a su b x c tetragonal unit cell reported in 
Ref . Q (both h and k are even) . Gray ticks show the super- 
structure reflections violating this unit cell (the h and/or k 
indices are odd) .— The arrow denotes the magnetic reflection 
emerging below Tjv — 7 K. 




FIG. 2. (Color online) Low-temperature crystal structure of 
(CuCl)LaTa2 07 with zigzag [CuCl] ribbons running along the 
b direction. The arrows show the displacements of oxygen 
atoms upon the a°b~c° tilting distortion (left panel) and the 
resulting displacements of Cu atoms (Aicu, right panel). The 
displacement Aa^cu is somewhat exaggerated for better visu- 
alization. The angle a measures the octahedral tilt, whereas 
ifi and £ characterize the Cu-Cl-Cu (Ji) and Cu-Cl-Cl-Cu 
(J4) superexchange pathways, respectively. Geometrical pa- 
rameters are listed in Tabic [Til 

even and odd h + k were observed (Fig. [TJ, hence the 
C-centered unit cell reported for the Br compound could 
not be applied to (CuCl)LaTa20 7 , and the primitive unit 
cell of (CuCl)LaNb20-;-i^— was used instead. Since the 
reflection conditions hOl, h = 2n and Okl, k = 2n are con- 
sistent with the Pbam space group of (CuCl)LaNb207, 
the respective crystal structure was introduced in the re- 
finement as the starting model. 

Despite the overall successful refinement, the atomic 
displacement parameter (ADP) of CI remained relatively 
high and possibly indicated the splitting of the CI posi- 
tion, as previously observed for CI atoms in single crys- 



TABLE I. Atomic positions and isotropic atomic displacement 
parameters Ui BO (in 10~ 2 A 2 ) for (CuCl)LaTa207 according 
to the refinement of the neutron data at 1.8 K. Lattice param- 
eters: a = 7.7663(5) A, b = 7.7640(3) A, c = 11.7374(5) A. 
Space group: Pbam. The Ui so of oxygen atoms were refined as 
a single parameter. The error bars are based on the Rietveld 
refinement. 





Position 


X 


y 


2 


C/iso 


Cu 


4ft 


0.7366(6) 


0.499(2) 


1 
2 


0.47(7) 


CI 


4/i 
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0.33(3) 



tals of (CuC^LaNbzOTjRef. [T3) and for Br atoms in 
(CuBr)LaNb 2 7 (Ref.GJ). However, the splitting of the 
CI position neither reduced the ADP nor improved the 
refinement. No signatures of out-of-plane displacements 
of the CI atoms were found either. 

Refined atomic positions and main interatomic dis- 
tances for (CuCl)LaTa207 are listed in Tables HI and [TTl 
respectively. Cu atoms have the typical four-fold coordi- 
nation (CUO2CI2 plaquette) with two short Cu-Cl bonds 
in the ab plane and two Cu-0 bonds along the c direc- 
tion. The short Cu-Cl bonds form zigzag ribbons along 
the b direction, whereas the TaOg octahedra develop a 
tilting distortion according to the a°b~c° pattern^ In 
(CuCl)LaTa207, all [CuCl] ribbons have the same orien- 
tation, whereas in (CuBr)LaNb20 7 two types of ribbons 
related by a mirror symmetry are disordered in the av- 
eraged crystal structure (the space group Cmmm) and 
form alternating layers along the c direction on the short- 
range scale (space group Ibam)^ A similar type of dis- 
order may be responsible for the somewhat high ADP of 
CI in (CuCl)LaTa207. However, the [CuCl] zigzag rib- 
bons of different orientation could be present as single 
defects, only. These defects do not manifest themselves 
as the split CI position in the averaged crystal structure 
probed by neutron diffraction. 

It is further instructive to compare details of the 
atomic arrangement in the (CuX)LaM207 scries. Ta- 
ble [TT] shows that the [LaM^C^] perovskite slabs fea- 
ture similar structures in all three compounds, with only 
a slight change in the tilting angle a. Main changes 
are observed in the [CuX] layers, where the larger size 
of Br causes longer Cu-X distances and, consequently, 
the smaller Cu-X-Cu angle ip compared to the CI com- 
pounds. There is also a slight, yet significant difference 
between the Cu-X-Cu angles in (CuCl)LaNb2C>7 and 
(CuCl)LaTa 2 7 3 This difference is confirmed by DFT- 
based structure relaxations^! that arrive at the Cu-Cl- 
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TABLE II. Selected geometrical parameters in the 
(CuX)LaM207 series: distances (in A) and angles (in deg). a 
denotes the tilting angle measured between the M-03 bond 
and the c direction, ip and £ describe the Cu-X-Cu (Ji) 
and Cu-X-X-Cu (J4) superexchange pathways, as shown in 
Fig. [2] Aa;cu (in A) measures the respective displacement 
of Cu atoms along the a axis with respect to x = \ and 
I (see Fig. [2]). Crystallographic data for (CuGl)LaNb 2 7 
and (CuBr)LaNb2C>7 are taken from Refs. [l7| (neutron 
refinement) andfisl. respectively. 

(CuCl)LaTa 2 7 (CuCl)LaNb 2 7 (CuBr)LaNb 2 7 
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FIG. 3. (Color online) High-resolution XRD patterns of 
(CuCl)LaTa 2 7 measured at 4.2 K and 750 K. The dashed 
line shows the Rietveld refinement of the 4.2 K data. Black 
ticks indicate the reflections of the a su b x a su b x c tetragonal 
unit cell reported in Ref. [l4| (both h and k are even). Gray 
ticks show the superstructure reflections violating this unit 
cell (the h and/or k indices are odd). Note that the supercell 
reflections are few and weak, and disappear at 750 K. The 
tiny peak remaining at 29 ~ 9.2 deg at 750 K is the subcell 
reflection 222. 



Cu angles of 109.6 deg and 108.4 deg in the Nb and Ta 
compounds, respectively. 

The slight variation in the Cu-Cl-Cu angle can be ex- 
plained by the different displacements of Cu atoms along 
the a axis (see Table [TTJ) . Previously, we have argue d 15 i 18 
that the Cu displacements are related to the tilting dis- 
tortion in the pcrovskite slabs, because the 03 atoms 
shift along the a direction upon the a°b~c° tilt and in- 
duce similar displacements of the Cu atoms, thus keeping 
the Cu-03 bonds perpendicular to the Cu-Cl bonds in 
the ab plane (see Fig. [2]). In (CuCl)LaTa207, the smaller 
tilt leads to the smaller displacement Axc u and, there- 
fore, results in the slightly reduced Cu-Cl-Cu angle. This 
mechanism elucidates the influence of non-magnetic pcr- 
ovskite slabs onto the magnetic [CuX] layers, and enables 
a deliberate modification of the magnetic layer via cation 
substitutions. 

To investigate temperature evolution of the 
(CuCl)LaTa20 7 structure, we performed high-resolution 
XRD experiments in the 4.2 — 750 K temperature range. 
In contrast to (CuCl)LaNb 2 7 and (CuBr)LaNb 2 7 , 
the observation of structural changes was impeded by 
the fact that the superstructure reflections are barely 
visible in the XRD data (Fig. |3|). Indeed, the super- 
structure originates from the displacements of light CI 
and O atoms, while heavy Ta atoms do not contribute 
to the respective reflections. Additionally, the weak 
orthoi'hombic splitting could not be resolved because 
of the reflection broadening. Since the independent 
evaluation of the a and b parameters led to a very 



unstable refinement, we constrained a = b in the whole 
temperature range. 

A comparison of the x-ray patterns measured at 4.2 K 
and 750 K demonstrates that the superstructure reflec- 
tions disappear upon heating. While the a lattice pa- 
rameter shows conventional thermal expansion with only 
slight changes in the slope, the peculiar temperature evo- 
lution of the c parameter signifies structural phase tran- 
sitions in (CuCl)LaTa20 7 (Fig. |4]). The increase in the 
c value below 450 K is followed by a sharp decrease be- 
tween 450 K and 650 K, and an eventual increase above 
650 K. These changes are likely unrelated to chemical 
transformations, because the sample weight is nearly un- 
changed up to 800 — 850 and the heating/cooling 
processes are fully reversible, thus ruling out the possi- 
bility of a decomposition. Unfortunately, the x-ray data 
are insufficient for a precise structure determination ow- 
ing to the very low intensity of the superstructure re- 
flections. High-temperature neutron experiments would 
be required to establish the structural changes associated 
with the transitions. Nevertheless, the app arent analogy 
to (CuCl)LaNb 2 07 (see Fig. 4 in Ref. suggests that 
the tilting distortion is eliminated upon the first struc- 
tural transformation at T\ ~ 450 K, while further heat- 
ing destroys the ordered arrangement of CI atoms above 
T2 — 650 K. These transition temperatures are similar to 
Ti ~ 500 K and T 2 = 620-640 K in (CuX)LaNb 2 7 ^^ 
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FIG. 4. (Color online) Temperature evolution of subcell lat- 
tice parameters of (CuCl)LaTa207. Different shadings denote 
the temperature ranges with dissimilar trends of the c param- 
eter. The transition at T2 — 650 K is also evidenced by the 
slight change in the slope of the a lattice parameter and of 
the subcell volume. Lines are guide-for-the-eye. 



IV. MAGNETISM 
A. Microscopic model 
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FIG. 5. (Color online) LDA density of states for 
(CuCl)LaTa207. The Fermi level is at zero energy. 




FIG. 6. (Color online) LDA band structure (thin light lines) 
and the fit with the tight- binding model (thick dark lines). 
The Fermi level is at zero energy. 



The close similarity between the crystal structures of 
(CuCl)LaNb 2 7 and (CuCl)LaTa 2 7 suggests that both 
compounds feature the same spin lattice with only a 
slight difference in the relevant microscopic parameters. 
To elucidate these subtle variations, we use DFT calcu- 
lations that evaluate individual exchange couplings. 

The LDA band structure of (CuCl)LaTa207 is shown 
in Figs. [5] and [6] The valence band comprises oxygen 2p 
states between — 7 eV and —1 eV as well as Cu 3d and CI 
3p states above —5 eV. The Fermi level is crossed by sev- 
eral bands forming a narrow complex between —0.3 eV 
and 0.3 eV. This band complex represents the d x 2_ y 2 
crystal- field levels of Cu +2 in the planar CUO2CI2 envi- 
ronment, where local x and y axes are directed along the 
Cu-0 and Cu-Cl bonds. Ta 5d bands are found above 
0.85 eV, compared to 0.4 — 0.5 eV for Nb 4d bands in 
(CuX)LaNb207. The spurious mctallicity of the LDA 
energy spectrum should be ascribed to the strong under- 
estimate of correlation effects in LDA. DFT+U calcula- 
tions result in realistic insulating spectra. 

Similar to Rcfs. [TH and[l9|, we use two complementary 
approaches to the evaluation of exchange couplings. The 
first approach is based on the tight-binding fit of the LDA 



band structure, with the resulting hoppings U introduced 
into an effective one-orbital Hubbard model featuring the 
on-site Coulomb repulsion U c g. Since ti <C U c s and the 
bands are half-filled, low-lying (magnetic) excitations can 
be described by a Heisenberg model with the purely AFM 
exchange J AFM = 4tf/U e $. The U cS = 4 eV value is used 
as a reasonable estimate for Cu +2 (see Refs. [33rl35T). 

The second approach involves local spin density ap- 
proximation (LSDA)+[/ calculations for a number of 
collinear spin configurations^ Resulting energies are 
mapped onto a classical Heisenberg model to obtain full 
exchange couplings J,. LSDA+[/ introduces a mean- 
field correction for correlation effects in the Cu 3c? 
shell, and provides accurate estimates of J, (see, e.g., 
Refs. [35|, H3, and HH) that are, however, rather sensitive 
to input parameters of the computational method. 

Previously, wc have argue d 33 ' 34 ' 39 that the on-site 
Coulomb repulsion parameter Ud depends on the double- 
counting-correction (DCC) scheme that subtracts the 
part of the correlation energy already contained in LSDA. 
For each of the DCC schemes - around-mean-field (AMF) 
and fully-localizcd-limit (FLL) - the Ud value has been 
adjusted with respect to the experimental Curie- Weiss 
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AMF solution FLL solution 



FIG. 7. (Color online) Left panel: magnetic couplings in the [CuCl] layer. Middle and right panels: the respective spin 
lattices according to the AMF and FLL solutions in Table ILTT1 In the middle panel, the shading shows the spin dimer formed 
in (CuCl)LaNb2C>7 by the coupling J4, while empty and filled circles denote the stripe AFM order observed experimentally 
in (CuCl)LaTa2C>7 and (CuBr)LaNb207. In the right panel, the shading indicates the Shastry-Sutherland plaquette with the 
intradimer coupling J4 and frustrated interdimer couplings J\ and J2, both FM. 



temperature and saturation field of (CuCl)LaNb207, to 
obtain U d = 4.5 eV and 8.5 eV in AMF and FLL, re- 
spectively (see Ref. [l^). Here, we use the same Ud val- 
ues, because the local environment of Cu +2 is essen- 
tially the same as in the Nb compound (see Table [TTJ) . 
In (CuBr)LaNb207 (Ref. fl8h . the adjusted Ud values are 
slightly different (5 eV and 12 eV, respectively), owing 
to the different nature of the ligands (Br instead of CI) 
and different coordination environment. The on-site ex- 
change parameter Jd is fixed at 1 eV. An elaborate dis- 
cussion on the role of the DCC in the evaluation of mag- 
netic couplings within LSDA+J7 can be found in Refs.l33l 
and|M 

Tabic Mil summarizes the computed exchange cou- 
plings for all three compounds of the (CuX)LaM2C>7 fam- 
ily. Here, we use J\, J{, and J" for inequivalent couplings 
between nearest neighbors, J 2 and J' 2 for couplings be- 
tween next-nearest neighbors, etc, and Jj_ for the inter- 
layer coupling along the c direction (Fig. [7]). The hop- 
ping parameters ti are obtained from fits to the LDA 
band structure using Wannier functions adapted to spe- 
cific orbital characters^ These hopping parameters are 
directly related to AFM exchanges J AFM , whereas the 
Ji values are an independent estimate of full exchange 
couplings that combine J AFM with ferromagnetic (FM) 
contributions Jf M . 

Basic features of the spin lattice of the (CuX)LaM207 
compounds have been extensively discussed in Refs. [l8| 
and Il9l Here, we only notice that the microscopic sce- 
nario is rather counter-intuitive because of the strongest 
AFM coupling between fourth neighbors. This unex- 
pected result was originally derived from the inelas- 
tic neutron scattering experiment on (CuCl)LaNb207 
(Ref. [Toh and later explained by the orthorhombic struc- 
tural model with the efficient, albeit long, Cu-CI-Cl-Cu 
superexchange pathway 16 i 19 i 42 Another salient feature of 
the (CuX)LaM2C>7 family is the sizable and nearly con- 
stant interlayer coupling J± mediated by the low-lying 
Nb 4c? and Ta 5d states. An experimental signature of 
this effect is the large hyperfine coupling for Nb atoms, 
as probed with nuclear magnetic resonance (NMR)4£ 



Despite certain discrepancies between the AMF and 
FLL results, exchange couplings evaluated by the same 
method enable a direct comparison to the structural data 
summarized in Table [II] The variation in the Cu-X-Cu 
angle ip has strong effect on the FM nearest-neighbor ex- 
change J\. The reduction in the angle from 109.0 deg 
in (CuCl)LaNb 2 7 to 108.2 deg in (CuCl)LaTa 2 7 and 
eventually to 101.8 deg in (CuBr)LaNb207 leads to a 
systematic increase in the absolute value of J\ , in agree- 
ment with Goodcnough-Kanamori rules. Note that both 

jAFM and jFM = j. _ jAFM arc c h an ging with the an _ 

glc. While J AFM depends on the orbital overlap, J™ is 
controlled by the Hund's coupling on the ligand site4i 

The leading coupling J 4 remains nearly unchanged 
upon the Nb/Ta substitution, in agreement with similar 
C angles pertaining to the curvature of the Cu-CI-Cl-Cu 
superexchange pathway (Table ILT1 and Fig. [5]). Although 
in (CuBr)LaNb207 the £ angle is notably decreased and 
the pathway becomes more curved, the intradimer ex- 
change J4 increases owing to the larger spatial extent of 
Br 4p orbitals compared to CI 3p. The AFM couplings 
J{, J", and J 2 are also enhanced. 

The AMF and FLL calculations arrive at similar mi- 
croscopic models that entail spin dimers formed by the 
leading coupling J 4 (Fig. [7J middle panel). The inter- 
dimer couplings include the FM nearest-neighbor interac- 
tion Ji, the AFM fourth- neighbor interaction J' A , and the 
AFM interlayer interaction J±. Next-nearest-neighbor 
couplings are also present, but their nature remains con- 
troversial. While AMF suggests a sizable AFM J2 and a 
weakly FM J 2 , FLL puts forward the sizable FM J 2 and 
a smaller J2 that may be either FM or AFM depending 
on the compound. The AMF- and FLL-based scenar- 
ios are, therefore, different with respect to the possible 
frustration of the spin lattice. The AFM J2 and negli- 
gible J'2 in the AMF-based model are compatible with 
other couplings that altogether establish the stripe or- 
dering pattern (parallel spins along 6, antiparallel spins 
along a and c, see the middle panel of Fig. [7J and the 
inset to Fig. [5]). The sizable FM J 2 in the FLL-based 
model would, in contrast, frustrate this magnetic order 
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TABLE III. Calculated exchange couplings in the (CuX)LaM2 07 series. For each compound, the first column lists hopping 
parameters of the tight-binding model U (in eV), the second column contains the derived AFM contributions to the exchange 
jAfm _ /jj eS ( m with U e g = 4 eV, while the third and fourth columns contain full exchange couplings Ji (in K) obtained 
from LSDA+(7 calculations with the AMF and FLL double-counting correction schemes. See text for details. The results for 
(CuBr)LaNb2C>7 are taken from Ref . [Isl The results for (CuCl)LaNb207 are calculated for the structural data from Ref . \vi\ 
and slightly differ from those reported in Ref. previously. 
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and make the situation more complex. 

Both AMF- and FLL-bascd models of (CuX)LaM 2 7 
have been considered in the literature. The AMF-based 
model (Fig. [7j middle panel) entails spin dimers with an 
intricate combination of non-frustrated interdimer cou- 
plings. Both thermodynamic and ground-state proper- 
ties of this model can be precisely evaluated by QMC, as 
shown for (CuCl)LaNb 2 7 (Ref. 01), (C uBr)L aNb 2 Q 7 
(Ref. [H), and (CuCl)LaTa 2 7 alike (Sec. llTCj) . 

The FLL-based model may be thought of as a modified 
Shastry-Sutherland spin lattice, where spin dimers are 
coupled by frustrated FM couplings J\ and J 2 forming 
triangles (Fig. [JJ right panel). This interpretation was 
proposed by Tassel et al.^ who performed FLL calcu- 
lations only and basically ignored the sizable AFM cou- 
plings J' A and J± that are both missing in the Shastry- 
Sutherland geometry. Despite a subsequent theoretical 
work^ the application of the Shastry-Sutherland model 
to (CuCl)LaNb 2 7 remains an obscure issue, because 
both thermodynamic and ground-state properties of the 
model, appropriately augmented with J4 and J±, arc 
hard to evaluate with a good precision. The QMC tech- 
niques fail, owing to the sign problem in a frustrated 
system, whereas exact diagonalization can not access suf- 
ficiently large finite lattices. Therefore, the application 
of the Shastry-Sutherland model to (CuCl)LaNb 2 7 is, 
at least presently, not amenable to the experimental ver- 
ification. 

In the following, we develop the microscopic descrip- 
tion along the lines of Refs. [IH and [HI, by focusing on 
the non-frustrated AMF-based model and exploring its 
properties with QMC. We further comment on the appli- 
cability of the Shastry-Sutherland model and its possible 
experimental verification (see Sec. IV)) . 




20 30 40 

Scattering angle, 20 (deg) 

FIG. 8. (Color online) Refinement of the magnetic structure. 
Experimental, calculated, and difference patterns are shown. 
Ticks denote the positions of magnetic reflections. The exper- 
imental and calculated patterns are offset from zero to ensure 
positive intensities. The inset shows the refined stripe AFM 
structure, with spins pointing along the b direction. 



B. Magnetic structure 

Kitada et aZ.— conjectured the stripe AFM ordering in 
(CuCl)LaTa 2 7 based on the experimental observation 
of the strongest magnetic reflection, which is labeled as 
10^ in Fig.[U Here, wc reconsider the magnetic structure 
of this compound in a powder neutron experiment that 
covers a wide angle range, thus giving access to all ob- 
servable magnetic reflections and improving the estimate 
of the ordered magnetic moment. 

The neutron diffraction pattern measured at 2 K re- 
veals one magnetic reflection only. The subtraction of the 
20 K pattern, where only nuclear scattering is present, 
provides a full picture of magnetic scattering, with few 
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additional weaker reflections at higher angles (Fig. [8]). 
This subtracted pattern is refined as a purely magnetic 
phase. All magnetic reflections could be indexed with 
the propagation vector k = (0, 0, \), which together with 
the Ah position of Cu allows for several irreducible rep- 
resentations. However, only one of these representations 
allowed the successful refinement. 

The resulting magnetic structure features stripes of 
parallel spins along either a or b direction, with antipar- 
allcl spin arrangement along the c direction and along the 
direction perpendicular to the stripes. Spins are directed 
along the stripes. Owing to the weak orthorhombic split- 
ting, the neutron data are not sufficient to decide whether 
the stripes run along a or b. However, the microscopic 
magnetic model gives clear indications for the stripes ar- 
ranged along the b direction according to the FM cou- 
pling Ji (Fig. [7]and Table [En]). This would also match 
the magnetic structure of (CuBr)LaNb207, where the ar- 
rangement of stripes along the [CuBr] zigzag ribbons (b 
direction) is confirmed in an NMR experiment ! 18 ' 46 

The refined magnetic moment fj, amounts 
to 0.53(1) fin, which is 0.2 [1b lower than in 
(CuBr)LaNb 2 7 {fJ, = 0.72(1) /i B , Ref. 0) and in 
the previous estimate by Kitada et al. (fi = 0.7(1) ^b)^ 
Phcnomenologically, the lower magnetic moment is well 
in line with the lower Neel temperature Tjv — 7 K 
compared to TV ~ 32 K in the Br compound (see 
also T/V/J4 in Table Hv|) . The lower magnetic moment 
and the reduced Neel temperature suggest that in 
(CuCl)LaTa207 quantum fluctuations are enhanced 
compared to its (CuBr)LaNb207 analog. 



C. Model simulations 

To quantify the magnetic model of (CuCl)LaTa207, 
we use QMC simulations for the Heisenberg spin Hamil- 
tonian and the spin lattice derived from the DFT cal- 
culations (Table ITTTl AMF solution). Unfortunately, the 
large number of inequivalent couplings prevents us from 
an independent evaluation of each exchange integral. We 
rather consider two relevant parameters, the intradimcr 
coupling J4 and the effective interdimer coupling J e g. 
The latter combines the interactions J\, J2, J4, and 
Jj_ , which are all of similar magnitude (see Table llllj) . 
Therefore, we assume — J± = J 2 = J 4, = J±_ = j and 
J c ff = 2| J\ \ + 2 J 2 + J4 + 2 J± = 7j, where the interdimer 
couplings are summed up according to their number per 
magnetic site (coordination number). 

The magnetic susceptibility of (CuCl)LaTa207 shows 
a broad maximum at T max ~ 10.5 K followed by the 
AFM ordering at Tjv — 7 K, which is hardly visible in 
the magnetization data (Fig. [9|) yet clearly identified by 
neutron measurements £ Since the susceptibility does not 
show exponential decrease at low temperatures, the spin- 
gap scenario established for (CuCl)LaNb 2 07 (Ref. \Wt ) 
can be safely excluded. The low-temperature behavior of 
the Ta compound is characteristic of a long-range ordered 
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FIG. 9. (Color online) Magnetic susceptibility of 
(CuCl)LaTa2 07 measured in the applied field of 0.5 T and 
the QMC fit with J4 ~ 16 K, J c g/Ji = 1.6. Deviations at 
low temperatures are likely related to an impurity contribu- 
tion. The inset shows the magnetization isotherm measured 
at 1.3 K (experimental data are taken from Ref. @) and the 
respective QMC fit for the same parameter set. 



TABLE IV. Magnetic properties of the (CuX)LaM 2 7 com- 
pounds: the Neel temperature (Tjv), ordered magnetic mo- 
ment fi, saturation field (noH s ), intradimer coupling J4, and 
the effective interdimer coupling J e ff (see text for details). 
The data for (CuCl}LaNb 2 07 and (CuBr)LaNb 2 7 are taken 
from Refs.Qjl and llH. respectively. 

(CuCl)LaTa 2 7 (CuCl)LaNb 2 7 (CuBr)LaNb 2 7 



Tjv/J 4 0.44 0.67 

fi (/Lts) 0.53 0.72 

fi H s (T) 22 30 84 

Ji (K) 16 25 48 

J cff /J 4 1.6 0.4 3.0 



quantum antiferromagnet that develops the short-range 
magnetic order evidenced by the susceptibility maximum 
at T max , and eventually forms the long-range ordered 
state at T N < T max . 

The microscopic magnetic model of (CuCl)LaTa 2 07 
should allow for the long-range AFM order at low tem- 
peratures. While the single susceptibility curve can be 
reproduced with different model parameters, the simul- 
taneous fit to the high-field magnetization isotherm taken 
from Ref. results in the unique solution with J4 ~ 16 K 
and Jes/Ji = 1.6 (Fig.[9]).i£ The fitted, powder -average 
g- value of g = 2.15 is in good agreement with the exper- 
imental estimates of g\\ = 2.13 and g± = 2.29 reported 
for (CuCl)LaNb 2 7I ii 

Similar to (CuCl)LaNb 2 7 and (CuBr)LaNb 2 7 , 
there is a certain ambiguity in the choice of interdimer 
couplings, because the overall coupling can be redis- 
tributed between weak AFM exchanges J2,Ji, and J±. 
Therefore, only the effective interdimer coupling is deter- 
mined with sufficient accuracy. The estimated J c s/ J4 ra- 
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FIG. 10. (Color online) Temperature dependence of the spin 
stiffness Lp s calculated for Lx Lx L finite lattices of different 
size L at J e s/Ji = 1.6. The shaded bar identifies the Neel 
temperature Tn/Ja — 0.34 corresponding to TV — 5.5 K for 
J 4 = 16 K. 



tios (Table lTV|) show clearly that the interdimer couplings 
are enhanced in (CuCl)LaTa207 compared to its Nb- 
containing counterpart, whereas (CuBr)LaNb207 fea- 
tures an even larger interdimer exchange. This evolu- 
tion of the spin lattice underlies the transition from the 
gapped spin-singlet ground state in (CuCl)LaNb 2 07 to 
the stripe AFM order in the other two compounds, as fur- 
ther discussed in Sec. [V] It also explains the increase in 
the Neel temperature and ordered magnetic moment in 
(CuBr)LaNb 2 7 compared to (CuCl)LaTa 2 7 (see Ta- 
bic HE}. 

The change in the intradimer coupling J4 is illustrated 
by the experimental values of the saturation field p-qH s 
(Table HVJ. Thus, the lower H s in (CuCl)LaTa 2 7 com- 
pared to its Nb-containing counterpart stems from the 
lower J4, even though the interdimer couplings are in- 
creased. The saturation field of (CuBr)LaNb 2 7 is, by 
contrast, about three times higher because of the larger 
J 4 and the pronounced increase in J e g. Altogether, our 
estimates of exchange couplings - J4 and J e g - are well 
in line with the experimental observations. 

To perform a more elaborate comparison between the 
proposed model and experiment, we estimated the or- 
dered magnetic moment p for the spin lattice with 
Jcg/ Ja = 1-6. Following the procedure described in 
Ref. Hll we computed static structure factors at the ex- 
perimental propagation vector k = (0, 0, i) for L x L x L 
finite lattices with L < 12, and performed the finite- 
size scaling. This way, the staggered magnetization 
m s = 0.681 \ib is obtained. To compare this model 
result with the experimental p, one should additionally 
take into account the spin-orbit coupling quantified by 
the g-value (fj, = gSm s with S = |) and the hybridiza- 
tion with the ligands. According to the LSDA+[7 calcu- 
lations, the hybridization reduces the magnetic moment 
by 25 - 30 %^ thus leading to \i = 0.51 - 0.55 p B in 
remarkable agreement with the experimental estimate of 
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FIG. 11. (Color online) Scaling procedure for locating the 
QCP. The spin stiffness (p„) multiplied by the squared lattice 
size (L 2 ) is plotted as a function of the effective interdimer 
coupling J e ff/J4 at the constant temperature of T/J4 = 0.02. 
The crossing point indicates the QCP at J c s/Ji ~ 1.14 sepa- 
rating the spin-singlet and long-range-ordered ground states. 
For the calculations, the L x L x L/2 finite lattices are used. 



0.53 [Ib- 

The Neel temperature can be evaluated from the tem- 
perature dependence of the spin stiffness p s . The spin 
stiffness is close to zero in the paramagnetic state (above 
Tjv) and reaches a finite value in the long-range-ordered 
state (below Tn). In the vicinity of Tn, a steep increase 
in p s is observed. The scaling properties^ of p s suggest 
that for a three-dimensional spin lattice the quantity Lp s 
is independent of the lattice size L at the transition tem- 
perature Tn- Therefore, Tn can be determined as the 
crossing point of Lp s (T) curves calculated for finite lat- 
tices with different L. The scaling procedure is shown in 
Fig. [TO]and results in T N /J A ~ 0.34 and T N ^ 5.5 K in 
good agreement with the experimental Tn — 7 K.— 



V. GENERAL MAGNETIC MODEL AND 
DISCUSSION 

The accurate microscopic description of 
(CuCl)LaNb 2 7 (Ref. [H ), (Cu Br)LaNb 2 Q 7 (Ref. M), 
and (CuCl)LaTa 2 7 (Sec. IIVC|) enables us to develop a 
common magnetic model for the (CuX)LaM 2 7 family. 
Although the three compounds slightly differ in the 
coupling regime (e.g., J[ and J" manifest themselves in 
the Br compound, while remaining weak in the Cl-bascd 
systems, see Table IIIII and Ref. [l8l ). their properties 
are well captured by the same spin lattice with the 
intradimer coupling J4 and the assumption of equal 
couplings — J 1 = J 2 = J'i = J±_ = j that contribute 
to the effective interdimer coupling J c ff = 7j. The 
magnetic behavior of this model is fully determined by 
the J e fi/ J4 ratio. At low J e ff, the spin dimers are weakly 
coupled and develop a spin gap in a spin-singlet (dimcr) 
state without the long-range magnetic order. Large 
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FIG. 12. (Color online) Magnetic phase diagram of the 
(CuX)LaM207 compounds. Triangles show calculated Neel 
temperatures Tjv for different J c g/J4, ratios and denote the 
boundary of the magnetically ordered phase. Dotted lines in- 
dicate the parameter regimes of three experimentally studied 
compounds, whereas large circles mark the experimental TV 
for (CuCl)LaTa207 and (CuBr)LaNb2C>7, where the stripe 
AFM order is observed. 



intcrdimer couplings close the spin gap and establish the 
stripe AFM order below Tjy. 

The regions of the spin-singlet ground state and stripe 
AFM order are separated by a QCP. The precise posi- 
tion of the critical point can be determined from simula- 
tions of the spin stiffness, similar to the estimate of the 
Neel temperature given in Sec. IIV CI The spin stiffness is 
close to zero in the spin singlet state, and reaches a finite 
value in the long-range-ordered state. Therefore, upon 
crossing the QCP the p s quantity changes abruptly, and 
for a three-dimensional spin system the crossing point of 
the L 2 p s (T) curves calculated for different lattice size L 
provides the position of the QCP at J e n/ Ji — 1-14 (see 
Fig. ED £1 

At low Jeff/Ji, the spin-singlet state continuously 
transforms into the high-temperature paramagnetic 
regime upon heating. By contrast, the elimination of 
the stripe AFM state is accompanied by an abrupt phase 
transition at T/v, which determines the phase boundary. 
The Neel temperature as a function of J e s/ Ji is eval- 
uated by QMC simulations for different J e s/ Ja ratios 
to obtain the phase diagram presented in Fig. [121 The 
(CuX)LaM 2 07 compounds are placed on this diagram 
according to their </ e ff/ J a ratios determined from QMC 
fits to the magnetic susceptibility and high-field magneti- 
zation (Fig. [9]). The comparison to the experimental T/v 
confirms the accuracy of our model that properly cap- 
tures the reduction in T/v in (CuCl)LaTa207 compared 
to (CuBr)LaNb 2 7 . 

Similar to TV, other parameters obtained from the ex- 
periment can be correlated with the J e g/Ji ratio. In 
Fig. [131 we plot the spin gap for the spin-singlet state at 
low J e fi/ J a and the staggered magnetization for the stripe 
AFM state at high J c s/ J±. Upon increasing the inter- 
dimer couplings, the spin gap is systematically reduced 
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FIG. 13. (Color online) Spin gap A (left axis) and staggered 
magnetization m s (right axis) in the microscopic magnetic 
model of the (CuX)LaNb207 compounds. Dotted lines indi- 
cate the parameter regimes of three experimentally studied 
compounds, whereas large circles are experimental results. 
See text for details. 



compared to A = J\ for an isolated dimer at J e s = 0, 
and eventually vanishes at the QCP. At higher J s/Ja, 
the staggered magnetization appears. It is worth noting 
that the trends for T/v and m s are somewhat dissimilar. 
While T/v depends on the overall energy of the exchange 
couplings and linearly increases with the interdimer cou- 
plings above J e ff/J4 ~ 1.5, the staggered magnetization 
approaches the saturated value of m s ~ 0.89 (is- This 
value is 11 % lower than the maximum magnetic moment 
of 1 hb for spin-i, and represents the effect of AFM spin 
fluctuations as well as residual quantum effects related 
to the low magnetic moment on the Cu site. 

Experimental spin gap and staggered magnetizations 
are in good agreement with our general model (see 
Fig. [T5|) . However, this comparison is less straight- 
forward than in the case of T/v, because m s is usually 
different from the ordered magnetic moment \i measured 
by neutron diffraction. Therefore, the experimental val- 
ues of /i should be scaled with the g- value and Cu-ligand 
hybridization, as explained in Sec. IIV CI While the ex- 
perimental fi value is not a direct measure of J e g/ J4, the 
changes in fi reflect the evolution of underlying exchange 
couplings and, e.g., manifest the reduction in T e ff/ Ji in 
(CuCl)LaTa 2 7 compared to (CuBr)LaNb 2 7 . The spin 
gap A can be estimated by different methods, with inelas- 
tic neutron scattering and high-field magnetization often 
leading to different results as, e.g., in (CuCl)LaNb2 07 J£ 
In our QMC simulations, we estimated A from the first 
critical field H c i, where the magnetization of the system 
starts increasing. Therefore, for the experimental esti- 
mate we used A ~ 14 K obtained from [i§H c \ ~ 10.5 T 
in (CuCl)LaNb 2 7 . 

Altogether, our general magnetic model elucidates the 
differences between the compounds of the (CuX)LaM 2 07 
family, and fully conforms to microscopic trends estab- 
lished by DFT calculations (Table ITTT)) . The replacement 
of Nb by Ta modifies the tilting angle of the MOg octahe- 
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dra within the [LaM^C^] slabs, thus slightly changing the 
Cu-Cl-Cu (ip) and Cu-Cl-Cl (£) angles in the ab plane, 
enhancing the FM coupling J%, and reducing the AFM 
coupling J4. The latter effect, although not well visible 
in the DFT results, is unequivocally established experi- 
mentally (Table HV|) . Although small at first glance, this 
alteration of the spin lattice is sufficient to change the 
magnetic ground state, because both (CuCl)LaNb207 
and (CuCl)LaTa207 are relatively close to the QCP in 
terms of the J c s/Ji ratio (see Figs. Q21 and fT5|) . Since 
(CuBr)LaNb2 07 lies further away from the QCP, the 
Nb/Ta substitution should have little effect on the mag- 
netism. Indeed, (CuBr)LaTa207 undergoes the AFM or- 
dering at Tn — 35 K, which is similar to TV — 32 K in 
(CuBr)LaNb 2 07.— 

Our results do not support the earlier conjecture by 
Kitada et at— who ascribed the qualitative difference be- 
tween (CuCl)LaNb 2 7 and (CuCl)LaTa 2 7 to a differ- 
ent regime of intcrlayer couplings mediated by the Nb 
and Ta atoms. We rather show that the interlayer cou- 
pling J± is essentially unchanged, whereas the different 
tilting angle of the NbC>6 and TaC>6 octahedra leads to 
slight changes in the [CuCl] layers, thus modifying the in- 
tralayer exchange couplings. In (CuX)LaM207, the vari- 
able magnetic behavior is related to structural changes 
in the ab plane, although the spin lattices are basically 
three-dimensional, with a sizable interlayer coupling Jj_ 
observed in all three compounds (Table HTT|) . 

The (CuX)LaM 2 07 compounds strongly resemble the 
family of spin-dimcr ACuX 3 halides (A = K, Tl; X 
= CI, Br), where quantum phase transitions between 
the gapped spin-singlet state and long-range AFM or- 
der were extensively studied experimentally. For ex- 
ample, TICuCbj features a spin gap A/ J ~ 0.6 at 
ambient pressure and transforms into the long-range- 
ordered antiferromagnet above 0.05 — 0.1 GPa.^ How- 
ever, in the ACuX 3 family the ordered ground state 
can not be stabilized at ambient pressure, unless an ad- 
ditional structural distortion is present, as experimen- 
tally found for NH4CUCI3 (Ref. (53|). An advantage of 
the (CuX)LaM 2 07 family is the possibility to stabi- 
lize both spin-singlet and long-range-ordered phases at 
ambient pressure. On the downside, single crystals of 
(CuX)LaM2 07 are exceedingly difficult to prepare^ and 
the exact pattern of intcrdimcr couplings has not been 
established experimentally. The spin lattice of ACUX3 
is similarly complex, yet amenable to the experimental 
study by inelastic neutron scattering on single crystals. 

The lack of a direct experimental information on 
the interdimer couplings resulted in a controversy re- 
garding the precise microscopic magnetic model of 
(CuCl)LaNb207 and the frustration of interdimer cou- 
plings in this compound (compare Ref s . [l6l and [ill) . Here, 
we elaborated on the non-frustrated version of the model 
(Fig. [3 middle panel),— and demonstrated its applicabil- 
ity to the whole (CuX)LaM2 07 family. This model not 
only captures the qualitative effect of different magnetic 
ground states depending on the interdimer couplings, but 



also provides a decent quantitative description, as shown 
in Figs. [T^] and Q21 Considering the enigmatic nature 
of early experimental results on (CuCl)LaNb207 and 
related compounds^ as well as futile attempts to un- 
derstand the underlying physics in terms of the square- 
lattice spin models^ our results are a remarkable ac- 
complishment that confirms the excellent potential of 
DFT and precise numerical simulations in the micro- 
scopic evaluation of complex spin systems. Neverthe- 
less, one has to keep in mind inevitable limitations of 
this approach that does not allow us to decide unam- 
biguously between the models with the frustrated and 
non-frustrated interdimer couplings (compare the right 
and middle panels of Fig. [7]). 

In Sec. IIV A[ we have argued that the FM Shastry- 
Sutherland model proposed by Tassel et alJ£ should be 
augmented by additional AFM interdimer couplings, and 
remains inaccessible to quantitative verification with fea- 
sible experimental techniques. Qualitatively, this model 
could be consistent with the quantum phase transition 
toward the stripe AFM ordering in (CuCl)LaTa207 and 
(CuBr)LaNb2 07, provided that one of the frustrating 
FM interdimer couplings becomes sufficiently weak or 
even AFM4£ Therefore, the magnetic frustration is not 
a generic feature of the whole (CuX)LaM2 07 family and 
may only affect the behavior of (CuCl)LaNb207, al- 
though even this conjecture requires further experimental 
verification. Future inelastic neutron scattering experi- 
ments on single crystals should elucidate to what extent 
the frustration and the Shastry-Sutherland-type physics 
are relevant for (CuCl)LaNb207. 

A natural reason for further experimental work on 
the (CuX)LaM2 07 family is the fact that the two 
Cl-containing compounds lie on different sides of the 
QCP. We have shown that a tiny structural effect re- 
lated to the Nb/Ta replacement is responsible for the 
change in the magnetic ground state. A similar change 
could be induced by hydrostatic pressure or chemical 
substitution. The available experimental data on the 
(CuCl)LaNb2- a; Ta :c 07 solid solutions^ suggest that the 
system separates into the long-range-ordered and spin- 
singlet phases. It would be interesting to understand 
whether or not this separation originates from a chemi- 
cal inhomogeneity, and whether or not an improved syn- 
thetic procedure could facilitate the continuous evolution 
of the system and the experimental access to the pro- 
posed QCP. 

In summary, we extensively characterized 
(CuCl)LaTa207 with respect to its crystal struc- 
ture, magnetic behavior, and microscopic magnetic 
model. Our results place (CuCl)LaTa207 between 
(CuCl)LaNb207 and (CuBr)LaNb207, and evidence 
sizable quantum fluctuations in this compound. The 
enhanced quantum fluctuations indicate the evolution 
toward the quantum critical point that separates the 
spin-singlet (dimcr) and long-range ordered ground 
states. The Nb/Ta replacement changes the tilting an- 
gles of the respective MOq octahedra and consequently 
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alters the positions of Cu atoms in the [CuCl] layers, 
thus enhancing the nearest-neighbor ferromagnetic 
coupling Ji, and reducing the intradimer coupling J4. 
This slight structural change is sufficient to induce a 
dramatic change in the magnetic ground state and drive 
the system across the quantum critical point. Therefore, 
both (CuCl)LaNb 2 7 and (CuCl)LaTa 2 7 as well as 
their solid solutions may be interesting materials for 
experimental studies of quantum critical behavior in 
spin-dimcr systems. 
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FIG. SI. Rietveld refinement of the D2B neutron data collected at 1.8 K. The excluded region 
around 29 = 40 deg is due to the cryostat window. Refinement residuals are Ri — 0.021, 
R p = 0.024, and R wp = 0.030. 
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FIG. S2. Thermogravimetric data for (CuCl)LaTa207. Note the onset of the weight loss 
around 850 K. 



